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ABSTRACT 

The binary J0737-3039A & B includes two pulsars in a highly relativistic orbit. The pulsed radio 
flux from pulsar B brightens considerably during two portions of each orbit. This phenomenon cannot 
be naturally triggered by the illumination of 7-rays or X-rays from pulsar A or the bow shock around 
pulsar B. Instead, we explain these periodic flares quantitatively as episodes during which pairs from 
pulsar A's wind flow into the open field line region of pulsar B and emit curvature radiation at radio 
frequencies within an altitude of '^ 10^ cm. The radio photons then travel through B's magnetosphere 
and eventually reach the observer on the other side of the pulsar. Our model requires that A's wind be 
anisotropic and that B's spin axis be somewhat misaligned relative to the orbital angular momentum. 
We estimate the expected 7-ray and X-ray emission from the system. 

Subject headings: pulsars: individual (PSR J0737-3039A, B) - stars: neutron - X-rays: stars 



1. INTRODUCTION 

In the double pulsar system PSR J0737-3039A & B dis- 
covered recently (Burgay et al. 2004; Lyne et al. 2004), the 
pulsed radio emission from B brightens by orders of mag- 
nitude near the orbital longitudes 210° and 280°, while 
during most of the orbit, B's radio emission is faint or 
non-detectable (Lyne et al. 2004; Ramachandran et al. 
2004). It has been heuristically suggested that B's ra- 
dio flares might be induced by the illumination of A, ei- 
ther in the form of particles, 7-rays, or radio photons 
(Lyne et al. 2004; Jenet & Ransom 2004). Assuming 
a distance of ^ 0.5 kpc (Lyne et al. 2004), the system 
was detected by Chandra to have an X-ray luminosity. 
La; - 2 X 10^° erg s"^ (McLaughlin et al. 2004). It is also 
within the EGRET point spread function of the unidenti- 
fied 7-ray source 3EG J0747-3412 (Hartmann et al. 1999), 
whose luminosity is '--^ 2 x 10'^^(Ai7-y) erg s^^ for the same 
distance and an unknown solid angle AfJ-y. 

2. BASIC EMISSION PROPERTIES 

We first estimate the emission properties of both pul- 
sars based on their measured period P and period deriva- 
tive P, by ignoring the interaction between them. Since 
the spindown power of A is much larger than that of B, 
Ea ^ 5.9 X 10^3 erg s^^ -> Eb ^ 1.6 x 10^" erg s^\ the 
interaction between the two pulsars is expected to have 
a negligible effect on A. The emission of B is, however, 
strongly influenced by the interaction, as we show in §3. 
• Pulsar A: The spin parameters Pa = 22.7 ms, and 
Pa = 1.74 X 10-1* (Burgay et al. 2004), yield a polar 
cap surface magnetic field Bp^A ~ 1-3 x 10^" G. The em- 
pirical relation between 7-ray luminosity L^ and E for 
known 7-ray pulsars (i.e. L^ oc E^/"^, Thompson 2003; 
Zhang & Harding 2000) yields L^^a ^ 1-3 x 10^3 erg s'^. 
More adequate pulsar acceleration models suggest a 7- 
ray luminosity in the range of (0.3-1) x lO^^ erg s~^ 
(Zhang & Harding 2000; Harding, Muslimov & Zhang 
2002). This estimate is consistent with the observed lu- 



minosity of the possibly related source 3EG J0747-3412, 
for Afi^ ~ (1 — 3). An associated X-ray luminosity 
Lx.a^^O^'^Ea ~ 5 X 10^° erg s-^ is expected from the 
full polar cap cascade (Zhang & Harding 2000), consis- 
tently with the empirical relation discovered in other spin- 
powered pulsars (Becker & Triimper 1997) and the Chan- 
dra data (McLaughlin et al. 2004). Like other millisecond 
pulsars, the radio emission of A has a wide double conal 
structure. Also the pulsar has a small inclination angle be- 
tween the spin and magnetic axes (Demorest et al. 2004). 
The Goldreich- Julian (1969) pair flux from both magnetic 

poles is Ngj,a = 2.8 x 10^2 s'^Bp^ioPI^ ~ 7 x lO^i s^^. 
• Pulsar B: The spin parameters are Pb = 2.77 s, and 
Pb = 0.88 X 10-15 (Lyne et al. 2004), yielding a polar 
cap surface magnetic field of Bps ^ 3.2 x lO^^ G. The 
pulsar polar gap is likely controlled by the IC process, and 
the 7-ray luminosity in this regime is L^^b '^ O.OSiEb '^ 
8 X 10^8 erg s'^ (Harding et al. 2002). The X-ray lumi- 
nosity may be estimated as ^ 10^^ Eb ^ 10^^ erg s^^. 
Both values are much smaller than the corresponding val- 
ues for A. The Goldreich- Julian flux from both poles is 
Ngj,b ~ 1.1 X 10^° s-i, and the expected pair multiplic- 
ity is Kb ~ OAkb,-i (Hibschman & Arons 2001). The 
intrinsic pair injection rate in the pulsar B magnetosphere 
is therefore 

N±{B) ~ 10^9 s-\ (1) 



3. PHYSICAL MECHANISMS FOR B'S RADIO FLARES 

Pulsar radio emission is broadly attributed to coherent 
emission by electron-positron pairs in the pulsar magneto- 
sphere. Our working hypothesis is that during B's radio 
flares, there is a signiflcant increase of the pair injection 
rate into the radio emission region. 

3.1. "f-ray and X-ray precipitations 

We consider it unlikely that B brightens because it is 
illuminated by A's radio beam, since radio waves do not 
lead to production of new e+e^ pairs. Jenet & Ransom 



(2004) suggested that the 7-ray beam may coincide with 
the radio beam of A and the energetic 7-ray photons would 
trigger a pair cascade in B's magnetosphere and lead to 
the radio flares. Putting aside the drawback that the 7- 
ray and radio beams are usually misaligned in known 7-ray 
pulsars (Thompson 2003), we demonstrate below that the 
physical parameters of this scenario are implausible. 

To maximize its effect, we assume that A's 7-ray lu- 
minosity is close to its highest expected value, L^^a '^ 
Xo33 g]-g g-i beaming into a solid angle Afi^ ~ 1. We 
also assume that the 7-ray spectrum resembles that of 
known 7-ray pulsars, i.e. the energy flux vF^ is flat be- 
tween 1 MeV and - 30 GeV (Thompson 2003). The range 
is bounded from above by the highest photon energy ca- 
pable of escaping from the polar cap cascade of A (see 
Eq. 29 in Zhang & Harding 2000). We may then express 
the injection spectrum as E'^N{E) — 10^^/log(3 x 10"') ^ 
10^^ erg s~^. In order to produce a pair in B's magneto- 
sphere, the energy of a 7-ray photon, E^, needs to satisfy 
{E-^/2meC^){B±/Bg) > x ^ 1/12 (Ruderman & Suther- 
land 1975), where B± = Bsin9kB and 9kB is the angle 
between the 7-ray momentum and the magnetic field B. 
Again, to maximize the effect we adopt sinOkB ^ 1- For a 
pure dipolc field, B ^ Bp{R/r)^ (where R is the stellar ra- 
dius and r is the radius at which pairs are produced), the 
threshold 7-ray energy is Ej^th ^ 1-2 McV(r/i?)^. Given 
E^^max, the maximum radius for pair production is then 
Tmax/R = (3 X lO^'/l. 2)1/3 ^ 29. The external pair cas- 
cade process is dominated by the synchrotron radiation 
of the higher generation pairs. Since the photon energy 
decreases by a factor 16 in each generation, each primary 
7-ray with energy E^ can generate 2^ pairs, where C, — 
[log(£^.^/£:.^,t,,)/log(16)] -h 1 (e.g. Zhang & Harding 2000). 
The radius-dependent solid angle (including both poles) 
of the open field line region is IS.Q.openi''') = An'^r/cP, and 
so the pair injection rate into B's open field region due to 
7-rays from A is 
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where the distance between the two pulsars is dAs ^ 
9 X 10^" cm (Lyne et al. 2004). This rate is negligible com- 
pared with the intrinsic pair injection rate of B itself (Eq. 
1). The total pair injection rate into B's magnetosphere 
(including the closed field line region) is 2 x 10^^ s"^ (calcu- 
lated by replacing Aflopen (?") by 47r in Eq. 2) , but it is be- 
lieved that pairs in the closed field line region can not con- 
tribute to the observed coherent emission from pulsars^. 

A's wind is terminated by the magnetic stress within 
B's magnetosphere through a bow shock. The distance 
of the bow shock from B is dsB = {Sfi'^d'AQc/ Ea)^^^ ^ 
5 X 10^ cm (Arons et al. 2004), where ^J,B ~ 3.75 x lO^'^ cgs 
is the magnetic moment of B. The bow shock produces a 
7-ray luminosity of ~ 10"^" erg s~^ at ^ 20 McV (Granot 



pair abundance in B's magnetosphere. The bow shock 
also produces X-rays, which may interact with the high 
energy 7-rays from the polar cap to produce pairs (see e.g. 
Zhang 2001; Harding ct al. 2002). However, for the esti- 
mated X-ray luminosity 10^^ erg s^^ (Granot & Meszaros 
2004), the number density of X-ray photons in the open 
field line region is only ^ 10^ cm^'^ (e^j/O.l keV)^^, lead- 
ing to an optical depth T^-y ~ 10^^ for a 7-ray traveling 
through the entire magnetosphere. The X-rays coming di- 
rectly from A have an effect that is smaller by a factor 

- (1030/1029) X [d^B/dABf - 0.03. 

We therefore conclude that the radio flares from B are 
not due to 7-ray or X-ray precipitations from either pulsar 
A or from the bow shock around pulsar B. 



3.2. Leakage of A's wind into B's magnetosphere 

The luminous A wind distorts B's magnetosphere into 
a shape analogous to the Earth magnetosphere as it is 
combed by the solar wind. According to the numerical 
simulation of Arons et al. (2004) , the vertical radius of the 
magnetospheric sheath is I ^ 7.5 x 10^ cm. Interpreting 
A's eclipse as synchrotron self-obsorption in the shocked 
A wind requires that the pair multiplicity of A be as high 
as KA -- 10^ (Arons et al. 2004; Lyutikov 2004). This 
multiplicity value is much larger than in standard cascade 
theory, <; 100 (Hibschman & Arons 2001). Hereafter, we 

normalize ka by 10®, although our principal conclusions 
remain valid at lower values. The total number of the pairs 
deposited to the bow shock region from A's wind is 



N±{A ^ sh) = Ngj.aka 



Afl-uj^A{dAB - dsB) 



^ 1.2 x 10^^ s-i KAA^^w^A/'iT:)-', (3) 

where Afl^A is the unknown solid angle of A's wind. 

It has been suggested that B spin axis aligns in the di- 
rection almost perpendicular to the orbital plane due to 
the external torque exerted by A's wind (Dcmorest et al. 
2004; Arons et al. 2004). Since the fine of sight (which 
sweeps across B's radio beam) is offset by 3° from the 
orbital plane (Lyne et al. 2004; Kaspi et al. 2004), B's 
magnetic axis must be oriented at a small angle relative 
to the orbital plane. It is therefore likely that A's wind 
would directly stream into at least one of the open fleld 
regions of B. The "leakage" may be realized through re- 
sistive effects (see Lyutikov 2004 and references therein). 
The analogy to the solar wind interaction with the Earth's 
magnetosphere suggests that the pairs gain access near 
B's magnetic pole. When this happens, some fraction of 
the pairs in A's wind can directly slide into B's mag- 
netosphere. At the same time the open field region in 
the "day" side facing A is greatly broadened due to the 
ram pressure of the wind. In principle, the stream from 
A's wind may encounter B's pair stream ("wind") along 
its path. The distance from B where the pressures of 
the two streams balance, dbB, can be found by equat- 



& Meszaros 2004). The ratio between the 7-ray flux from ing riBEB/{cAVLw.Bd\ 



bB) 



EAl[cArtw^A{dAB - ddBf 



A and this component is ^ (lO^^/lO^") x {dsB/dABY ~ 3, 
and so this component does not increase significantly the 

^When estimating the pair injection rates (Eqs. 1, 2), a pure dipolar geometry was assumed for B. In reality, B's magnetosphere is severly 
distorted by A's wind (Arons et al. 2004; Lyutikov 2004), but our conclusion remains unchanged when the distortion effect is included. 



where AVl,ui.B oc rffcs is the solid angle of B's wind, and 
B's wind luminosity is smaller than the spindown power 



by a factor of 773 < 1 due to pair screening. We get 

dbB -^ 1-5 X 10^ cm r]]^^{An^^A/^^w,By^^ < 4b- Once 
pairs from A's bow shock leak into B's magnetosphere, 
they will stream all the way down to B's surface due to 
the external pressure gradient from above. 

We parameterize the pair injection rate into B's open 
field line region as 

N±{sh ^ B) - 1.2 X lO^'^ s"i ?7_iKA ei^^w A/47r)"\ 

where rj = 0.1?7_i is the fraction of all pairs from A's bow 
shock that enter this region. The resulting injection rate is 
orders of magnitude larger than the intrinsic pair injection 
rate from B in equation (1), even if ka,6 is smaller than 
unity. We therefore suggest that these pairs are the cat- 
alyst for B 's radio flares. A natural coherent mechanism 
would be the two-stream instability between the down- 
stream A wind and the upstream B wind. The radio pho- 
tons emitted by these pairs travel downwards, but owing 
to the curved dipole geometry of the magnetic field that is 
responsible for their production, they are not blocked by 
the neutron star. Rather, the radio waves pass through 
the inner magnetosphere and the gravitational field of B 
and eventually reach the observer on the other side of the 
neutron star. 

The observed double conal emission of B is separated 
by about 0.03 s (Lyne et al. 2004). Since the open field 
region is nearly oriented towards the observer, the radio 
beam opening angle is 6b ^ 27r(0.03/2.77) ^ 0.07. Based 
on the dipole field geometry and the gravitational bend- 
ing effect, this angle should be 9b ~ 2{9t — Og), where 
6»j = (3/2)6'r = (3/2) (r/rL)i/2 is the angle between the 
tangent of the field line at the emission point and the 
magnetic polar direction, and 6g = 2rsch/b is the gravita- 
tional deflection angle. Here rsch = 2GM/c^ ^ 3.7 x 10^ 
cm is the Schwarzschild radius of the neutron star for a 
mass of 1.25 M© (Lyne et al. 2004), and h ~ riOt - Or) 
is the photon's impact parameter from the star. Due to 
A's wind pressure, the field lines on the "daytime" side 
are broadened. We still assume a dipolar geometry, but 
take the light cylinder radius r^ to be comparable to the 
vertical size of the magnetic sheath, r^ ^ I ^ 7.5 x 10^ 



1/2 



cm. This gives 9r ^ 0.1 2r, 
0.13rg , where rg — r/10^ 



0.17ro and 



2(0.17rg/^ - 0.13,8 



-3/2x 



We therefore obtain 

or Ob r^ 0.08 for rg - 1. 
Thus rg ~ 1 is required from the data for such a geometry. 
The average Lorentz factor of the pairs that enter 
B's magnetosphere could be estimated from the relation 
j±meC^ Ngj,aka = Ea- In pulsar magnetospheres, the 
synchrotron cooling time scale for pairs is very short so 
that the perpendicular component of the pair energy is 
lost instantaneously. The parallel Lorentz factor of the 
pairs is 7-1- y — ^||7± ^ 210k^ g^y, where ^|| < 1 is a geo- 
metric factor depending on the incidence angle of the pairs 
relative to the magnetic field (Zhang & Harding 2000). At 
a height r ^ lO^rg cm, the curvature radius of the field 
line for r^ ^ Z is p ~ 1.2 x lO^rg cm, and the typical 
curvature radiation frequency of the pairs is 



2^ 



3 c3^-l/2 



350 MHz K^%qr^ 



(5) 



The characteristic plasma frequency for the two stream 
instability is cUp = ^^^AATrn'e^ /m^Y''^ (Ruderman & 
Sutherland 1975; Medvedev & Loeb 1999), where n' - 
N±{sh — > i?)/7-|- ||7rr^0^c is the co-moving down-stream 



plasma density. If 



< 



0.25C 



5/6 l/3>, 

II 



(requiring ka,i 



< 



then the two-stream instability greatly am- 



plifies the curvature radiation (Ruderman & Sutherland 
1975). Given the uncertainties in the value of ka,6 a-nd 
^11 , the typical enhanced curvature radiation frequency is 
consistent with the 680-3030 MHz band at which B is de- 
tected (Lyne et al. 2004). We therefore conclude that the 
emission altitude is ~ 10* cm, a value consistent with the 
emission altitudes derived for other pulsars (Kijak & Gil 
2003). We note that although some models of pulsar ra- 
dio emission (e.g. Lyutikov, Blandford & Machabeli 1999) 
disfavor emission from regions near the magnetic pole, the 
observed narrow beam calls for this geometry, and we spec- 
ulate that the two-stream instability discussed above plays 
an essential role in achieving the strong coherence of the 
emission. 

3.3. Consequences of the model 

There are several direct consequences of our model. 
First, the pairs eventually deposit energy onto B's polar 
cap region, making B a strong X-ray emitter as well. The 
pairs lose energy via curvature radiation and IC as they 
stream towards B's surface, and both braking processes 
become more significant near the star surface. The cur- 
vature radiation power Jcr^tT'cC^ ~ ^(2c/3)(e^//9)7^ r^ 



10 



-7 



erg s 



=-1. 



-1/2 



is negligible. The IC braking is, how- 
ever, important. With 7-1- ^ 210, the IC process is not ini- 
tially in the resonant regime, i.e. ^±ex > {B.p/ Bq)mec'^ ~ 
35 keV(i3p/3 x 10^^ G), where ex is the typical X-ray pho- 
ton energy (Eq. 7), and Bq = 4.414 x 10^^ G is the critical 
magnetic field. However, even in the non-resonant regime, 
the IC power jcRfneC^ = — (4/3)7^ccrT^p/i is large for 
typical parameter values, so that the pairs undergo signif- 
icant deceleration at r < 3R. The IC process then enters 
the resonant regime, and the braking power becomes even 
larger. As the pairs hit B's surface, their typical Lorentz 
factor is of order a few. However, the IC 7-rays are beamed 
towards the surface, so that their energy is deposited on 
the near polar cap region as well. We therefore assume 
that all the initial pair energy is deposited on the polar 
cap region and later radiated as X-rays. The X-ray lumi- 
nosity can then be estimated as 

Lx,B '^ N±{sh -^ i?)7±j|meC^ 

^ 2.1 X 10^" erg s^^ ri^i^iAn^^A/^Tr)-^ , (6) 

yielding a value that is consistent with the Chandra ob- 
servations (McLaughlin et al. 2004) for 77 ^ 0.1 during 
the phases when B is illuminated by A's wind. The typi- 
cal temperature of the polar cap is Tpc ^ (Ll/aTrR^Y'^ '-^ 
3.1 X 10*^ Kr/i/;*Cy^'*(Afi^,A/47r)-i/4, with a typical photon 
energy 



2.8fcr, 



1/4^1/4. 



'pc ^ 0.7 keV 77^7^11 (Ar!^,A/4^) 



\-l/4 



(7) 



This X-ray component should be modulated at B's spin 
period of 2.77 s, but may not be modulated with orbital 



phase due to the long thermal inertial time of neutron stars 
(Eichler & Cheng 1989). 

At the higher altitudes of r ~ lO^cm, the IC pho- 
tons produced by the pairs will not be blocked by the 
neutron star and can reach the observer. Their typi- 
cal energy is e^ ~ l±^x ~ 30 MeV. The number den- 
sity of X-ray photons is rix ^ Lx^pd {^x'^t^t'^ c) ^ 2 x 
10^^ cn\~^riJ_^ ry , giving a Thomson optical depth tjc — 
nxTax ^ 1.3 X \Q^^r]\ Vj . The 7-ray luminosity is 

L^^B' ^ 30 MeV X N{sh -> B) x nc - 7.5 x 10^^ erg s-\ 
much lower than the intrinsic 7-ray luminosity of B, which 
itself is much lower than the 7-ray luminosity of A. 

Since we have attributed B's radio flares to the inter- 
action between A's wind and the open field region of B, 
the geometric model of Jenet & Ransom (2004) which is 
based on A's radio beam illuminating B, has to be modi- 
fied (keeping in mind that A's wind should be wider than 
its radio beam, allowing for a temporal mismatch between 
the two as reported [Ransom et al. 2004]). In order to ob- 
tain two specific orbital phases during which B flares, we 
require that A's wind be anisotropic, and that B's spin 
axis be somewhat misaligned relative to the orbital an- 
gular momentum. The lack of a detection of inter-pulses 
supports the latter assumption. The suggestion of two ad- 
ditional orbital episodes of enhanced weak emission (Ra- 
machandran et al. 2004) can be interpreted as the inter- 
action of A's wind with B's second magnetic pole (which 
has the less favorable orientation). For the two bright 
phases, the roughly-equal double-component radio profile 
near longitude 280° (Lync et al. 2004) is interpreted as A's 
wind streaming directly into B's open field region, while 
the weak precursor followed by an intense main pulse at 
longitude 210° (Lyne et al. 2004) can be explained if A's 
wind only covers partially B's open field region. For the 
other two dim phases (Ramachandran et al. 2004), either 
the wind interaction region only covers a small fraction of 
the open field region, or the line of sight only grazes the 
emission region. A more refined geometric model could be 
developed from additional radio data, including the con- 
sequences of geodetic precession (Jenet & Ransom 2004). 

In our model, the radio waves travel through near- 
surface closed field regions where the magnetic fields 
are stronger and the plasma is denser. At closest ap- 
proach to the neutron star, the local magnetic field is 
^ 1.6 X 10^" G. The local plasma frequency is utp^c ~ 



(47rn±^ce^/7±,cWe)^/^ ^ 110 MHz k^ ^j^ ^21 where n 



'±,c, 



kb,c — n±^c/ncj, and 7±^c are the number density, multi- 



plicity and typical Lorentz factor of the pairs in the closed 
field region (with these pairs being seeded by the 7-rays 
emanating from A, the bow shock, and the radio emission 
region) . The plasma cutoff frequency is much lower than 
the observed frequency, but could introduce a large disper- 
sion measure for B. The propagation through the plasma 
may also lead to novel polarization signatures. 

4. SUMMARY 

We have interpreted the periodical radio re-brightening 
of B as episodes during which pairs from A's wind flow 
into the open field line region of B and emit curvature ra- 
diation at an altitude of ^ 10® cm. The pair multiplicity 
required to explain A's eclipse as synchrotron absorption 
in the bow shock around B's magnetosphere, yields the re- 
quired radio frequency for the curvature radiation in B's 
magnetoshpere (Eq. 5) . The radio photons travel through 
B's magnetosphere and gravitational potential and even- 
tually reach the observer on the other side of the pulsar. 
Our model requires that A's wind be anisotropic, and that 
B's spin axis be somewhat misaligned with the orbital an- 
gular momentum. 

The system has several components of high energy 
emission. The dominant 7-ray source is A, whose lu- 
minosity Lj^A ~ (0.3 — 1) X 10^"^ erg s~^ is consistent 
with that of' the unidentified EGRET source 3EG J0747- 
3412. The source is expected to have a 22.7 ms period 
which could be detected by the future GLAST mission 
(http://glast.gsfc.nasa.gov/). In the X-ray band, there 
are three sources that could account for the luminosity 
measured by the Chandra satellite, L^ ^ 2 x 10'^" erg s^^. 
They are: (i) polar cap heating and cascade emission from 
A (Zhang & Harding 2000); (ii) B's polar cap heating by 
A's wind (this work); and (Hi) emission by the interstellar 
medium shocked by A's wind (Granot & Meszaros 2004). 
The first two components should be pulsed with the corre- 
sponding pulsar periods (22.7 ms and 2.77 s), respectively, 
and should possess a thermal spectral component. Fur- 
ther observations with longer exposure times are needed 
to separate these components. 
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